Evidence for the Direct Two-Photon Transition from i/){3686) to J/i/) 
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The two-photon transition ■(/'(SGSG) — >■ 77J/1/J is studied in a sample of 106 miUion ?/;(3686) 
decays collected by the BESIII detector. The branching fraction is measured to be (3.1 ± 
0.6(stat)lj g(syst)) x 10"** using J/ip — >■ e^e~ and J/ip — > fJ-'^ iJ-~ decays, and its upper limit is 
estimated to be 4.5 x 10~* at the 90% conference level. This work represents the first measurement 
of a two-photon transition among charmonium states. The orientation of the ^/'(3686) decay plane 
and the J/^ polarization in this decay are also studied. In addition, the product branching fractions 
of sequential El transitions 1/1(3686) "fXcJ, Xcj — > "^J /ijj^J — 0, 1, 2) are reported. 
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The XYZ [l| particles, which do not fit potential model 
expectations in QCD theory, have been a key challenge 
to the QCD description of charmonium- like states [3]. 
To fully understand those states, it is necessary to con- 
sider the coupling of a charmonium state to a DD me- 
son pair. These coupled-channel effects, which also play 
an important role in the charmonium transitions of low 
lying states {i.e., from -0(3686) to J/ip), are especial- 
ly relevant for the radiative transition processes 3]. In 
the well-known electric dipole transitions, the strength 
of coupled-channel effects will likely be hard to estab- 
lish, since the accompanying relativistic corrections may 
be more important 14!]. However, the two-photon transi- 
tion -0(3686) — > 'Y^J/ip is more sensitive to the coupled- 
channel effect and thus provides a unique opportunity to 
investigate these issues 

Two-photon spectroscopy has been a very powerful 
tool for the study of the excitation spectra of a variety 



of systems with a wide range of sizes, such as molecules, 
atomic hydrogen and positronium Q . Studying the anal- 
ogous process in quarkonium states is a natural exten- 
sion of this work, in order to gain insight into non- 
perturbative QCD phenomena. But so far, two-photon 
transitions in quarkonia have eluded experimental ob- 
servation 0-ll|- For example, in a study of 0(3686) — >■ 
lXcj{J = 0, 1, 2) reported by CLEO-c Q, the upper limit 
for S(0(3686) — )■ 77^7-0) was estimated to be 1 x 10"'^. 

This Letter presents the first evidence for the two- 
photon transition -0(3686) — ?► 77J/0, as well as studies 
of the orientation of the 0(3686) decay plane and the 
J/ip polarization in the decay. The branching fractions 
of double El transitions ?/'(3686) l{lJ l'4')xc.j through 
Xcj intermediate states are also reported. The data ana- 
lyzed were obtained by the BESIII experiment ^IQ] view- 
ing electron-positron collisions at the BEPCII collider. 
An integrated luminosity of 156.4 pb~^ was obtained at 



3 




0.2 0.4 0.6 0.2 0.4 0.6 0.2 0.4 0.6 



MT,(GeV/c') 

FIG. 1. Up (a): distributions of M^-^-rccoii in data (points) 
and in the combined dataset (solid line) of MC simulation 
of i/j(3686) decays (shaded histogram) and continuum back- 
grounds (dashed line), before the KF is applied. The arrows 
indicate the window to select a J/i/; candidate; Down: scatter 
plots of M^gjjj_recoii vcrsus M-y^ for the 776+6" channel, in 
data (b), continuum data (c), MC simulated signal (d), after 
applying the KF constrain and the A/^-y_recoii window. The 
corresponding plots for the 77/^+/^" channel are very similar. 

a center-of-mass energy ^/s = M('0(3686)) = 3.686 GeV. 
The number of '0(3686) decays in this sample is esti- 
mated to be (1.06 ± 0.04) x 10^ In addition, 42.6 
pb~^ of continuum data were taken below the "0(3686), 
at y/s — 3.65 GeV, to evaluate the potential backgrounds 
from non-resonant events. 

The upgraded BEPCII [12] at Beijing is a two-ring 
electron-positron colhder. The BESIII detector ^Oi] is 
an approximately cylindrically symmetric detector which 
covers 93% of the solid angle around the collision point. 
In order of increasing distance from the interaction point, 
the sub-detectors include a 43-layer main wire drift cham- 
ber (MDC), a time-of-flight (TOP) system with two lay- 
ers in the barrel region and one layer for each end- 
cap, and a 6240 cell CsI(Tl) crystal electro-magnetic 
calorimeter (EMC) with both barrel and endcap sections. 
The barrel components reside within a superconducting 
solenoid magnet providing a LOT magnetic field aligned 
with the beam axis. Finally, there is a muon chamber 
consisting of nine layers of resistive plate chambers with- 
in the return yoke of the magnet. The momentum resolu- 
tion for charged tracks in the MDC is 0.5% for transverse 
momenta of 1 GeV/c. The energy resolution for showers 
in the EMC is 2.5% for 1 GeV photons. 

This work studies ■0(3686) ^^Jjil) followed by 
J/'0 — > {£ denotes e or /i), which is referred to 

as the signal process. Events selected contain exactly 
two oppositely charged good tracks in the MDC track- 
ing system, corresponding to the dilepton from J/tp de- 
cay. The requirements to judge a track as good include 
I cos 6'! < 0.93 {9 is the polar angle with respect to the 
beam direction), and the minimum distance of approach 
between the track and the production vertex less than 
10 cm along the beam axis and less than 1 cm projected 
in the perpendicular plane. The lepton is identified with 
the ratio of EMC shower energy to MDC track momen- 



tum, E/p, which must be larger than 0.7 for an electron, 
or smaller than 0.6 for a muon. To suppress non-J/ip 
decay leptons, we require the momentum of each lepton 
to be larger than 0.8 GeV/c. A vertex fit (VF) constrains 
the production vertex, which is updated run-by-run, and 
the tracks of the dilepton candidates to a common vertex; 
only events with Xyp/'^-O-f- < 20 are accepted. 

Reconstructed EMC showers unmatched to either 
charged track and with an energy larger than 25 MeV in 
the barrel region (|cos0| < 0.80) or larger than 50 MeV 
in the end-caps (0.86 < | cos9\ < 0.92) are used as photon 
candidates. To reject bremsstrahlung photons, showers 
matching the initial momentum of either lepton within 
10° are also discarded. Showers from noise, not originat- 
ing from the beam collision, are suppressed by requiring 
the EMC cluster time to lie within a 700 ns window near 
the event start time. 

Events are required to have only two photon candi- 
dates. A kinematic fit (KF) constrains the vertexed 
dilepton to the nominal mass of the intermediate J/ip, 
and the resulting J/^p and photon candidates to the 
known initial four-momentum of the ■0(3686). The KF fit 
quality Xkf required to be x|-p/d.o.f. < 12. For con- 
venience, we use 7ig (7sm) to denote the larger (small- 
er) energy photon. As indicated in Fig. [TJa), J/ip 
candidates are identified with the requirement that the 
recoil mass of the two photons, Mjj^^ocoih is within 
(3.08,3.14) GeV/c2. 

Scatter plots of recoiling mass M^^^-rocoii from the 
lower energy photon 7si„ versus invariant mass of two 
photons Mj^ are shown in Fig. [1] where clear resonance 
bands are seen from the decays -0(3686) — )■ JXcjiJ = 
0, 1, 2) (three horizontal bands) and V'(3686) — > 7r°(?7) J/'0 
(two vertical bands). As indicated in Fig. [IJc), the con- 
tinuum backgrounds are most dominant at the top of the 
plots, of which the primary sources include the bhabha 
scattering, the dimuon process and the ISR production 
of J/-0. These backgrounds are excluded by discarding 
events with M-y^_^ -j-ccoii > 3.6 GeV/c^. To suppress back- 
grounds from V'(3686) — n'^ {rj) J / -ip , the diphoton invari- 
ant mass is required to be larger than 0.15GeV/c^ 
and the recoil momentum of the diphoton must be larger 
than 0.25GeV/c. 

Monte Carlo (MC) simulations of "0(3686) decays are 
used to understand the backgrounds and also to estimate 
the detection efficiency. At BESIII, the simulation in- 
cludes the beam energy spread and treats the initial-state 
radiation with KKMC [lj|. Specific decay modes from 
the PDG [145 are modeled with EVTGEN and the 
unknown decay modes with Lundcharm [16J. The detec- 
tor response is described using GEANT4 [l7|. For the 
■0(3686) — >■ jjJ/ip channel, the momenta of decay par- 
ticles are simulated according to the measured polariza- 
tion structure in this work. Generic ^/'(3686) decay sam- 
ples serve for understanding the background channels; 
dominant backgrounds were generated with high statis- 
tics. Angular distributions of the cascade El transitions 
■0(3686) — )■ 7Xcj —>■ 77 •//'0 are assumed to follow the 
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FIG. 2. (color online) Plot a: unbinned maximum likelihood 
fit to the distribution of M.^^^ -recoil in data with combina- 
tion of the two J/tp-deca,y modes. Thick lines are the sum of 
the fitting models and long-dashed lines are the Xcj shapes. 
Short-dashed lines represent the two-photon signal process- 
es. Shaded histograms are i/)(3686)-decay backgrounds (yel- 
low) and non-?/)(3686) backgrounds (green), with the fixed 
amplitude and shape taken from MC simulation and contin- 
uum data. Plot b: the number of standard deviations, n„, 
of data points from the fitted curves in plot a. The rates of 
the signal process and sequential Xcj processes are derived 
from these fits. Plot c: distributions of M-,-y_rocoii in data 
(signals and known backgrounds) with the kinematic require- 
ment 3.44GeV/c^ < M.,^^^-rccoa < 3.48GeV/c^ and with the 
removal of Xkf s^nd M.^.y_rccoii restrictions. Plot d: stacked 
histograms of the three Xcj components in plot c. 

formulae in Ref. [l^ ■ Note that the Xcj line shapes were 
simulated with the Breit-Wigner distributions weighted 
with E^tE^t to account for the double El transitions, 

'1 '2 

and extended out to ±200 MeV/c^ away from the nomi- 
nal masses, using masses and widths in PDG Here, 
E^* (£'7^) is the energy of the radiative photon 7^ (72) in 
the rest frame of the mother particle V'(3686)(xcj)- 

The yield of the signal process -0(3686) — > 777/0, 
together with those of the cascade El transition pro- 
cesses, is estimated by a global fit to the spectrum of 
M.y^^_iocoii- The fit results are shown in Fig. [2] The 
shape and magnitude of '0(3686)-decay backgrounds were 
fixed based on MC simulation. Non-'0(3686) decay back- 
grounds are estimated in continuum data, scaling by lu- 
minosity and the 1/s dependence of the cross sections. 
This scaling is verified by the good description of the 
J /if) backgrounds in the M^^-rocoii distribution shown in 
Fig. [TJa). The distributions of the signal process and 
the cascade El process are taken from the reconstructed 
shapes in MC simulation of the modes and smeared with 
an asymmetric Gaussian with free parameters, which is 
used to compensate for the difference in line shape be- 
tween MC and data. By taking the MC shape, detector 
resolution and wrong assignment of the El photon are 
taken into account. The quality of goodness-of-fit test, 
X^/d.o.f.= 108.0/94 = 1.15 in the 776+6" mode and 
124.8/94 = 1.33 in the 77//+//" mode. The observed sig- 
nal yields are given in Table HI The V'(3686) 11 J 1^ 
transition is observed with a statistical significance of 
6.6(7, as determined by the ratio of the maximum likeli- 
hood value and the likelihood value for a fit with nuU- 



TABLE I. For different channels: the number of observed sig- 
nals rte (rip) and detection efficiency ee (e^) in 776+6" (77/^+^") 
mode; the absolute branching fractions. On the bottom, the rel- 
ative branching fractions Rmn = ^Xcm /^Xcn ' where B^^j = 
B{rp{3686) —>■ 'y{'yJ/'4')xc i) listed. Here the first errors are 
statistical and the second are systematic. 



Channels nl 7J%) ^^A%) B{xW^ 

77J/i/' 564±116 22.4 536±128 30.0 3.1 ± O-el^-Q 

7(7J/V)xco 1801±60 19.3 2491±69 26.0 15.1 ± 0.3 ± 1.0 

jhJ/i>)x^i 59953±253 28.5 81922±295 38.2 337.7 ± 0.9 ± 18.3 

7(7J/j/i)x,, 32171±187 27.5 44136±219 37.1 187.4 ± 0.7 ± 10.2 

R21 = ^{%) Roi ^ ^(%) R02 = |^(%) 

55.47 ± 0.26 ± 0.11 4.45 ± 0.09 ± 0.18 8.03 ± 0.17 ± 0.33 



signal hypothesis. When the systematic uncertainties are 
taken into account with the assumption of Gaussian dis- 
tributions, the significance is evaluated to be 3.8(7, which 
corresponds to a probability of a background fiuctuation 
to the observed signal yield of7.2xl0"^. The upper lim- 
it for S(0(3686) 77 J/V-) is estimated to be 4.5 x lO"-* 
at the 90% confidence level including systematic uncer- 
tainties. 

In calculating BdjJ/ip), a correction factor is includ- 
ed due to the interferences among Xcj states. This ef- 
fect was checked by the variations of the observed signals 
in the global fit with inclusion of a floating interference 
component, which is modeled by the detector-smeared 
shape of a theoretical calculation [5]. It is found that 
relative changes on the signal yields are negative with 
lower bound of —10%. Hence, a correction factor 0.95 is 
assigned and 5% is taken as systematic uncertainty. 

A cross-check on our procedures is performed with 
the M-y.y_recoii spcctrum for the events in the region 
3.44GeV/c2 < M^^.„_rccoii < 3.48GeV/e2 without re- 
strictions on Xkf ^^"^ A/.77-recoii, as shown in Fig. [2Jc). 
An excess of data above known backgrounds can be seen 
around the J/ip nominal mass, which is expected from 
the sought-after two-photon process. With the inclusion 
of the estimated yields of the signal process, the excess 
is well understood. The high- mass peak above the J/V^ 
peak comes from the backgrounds of ?/'(3686) — 'k^'k^J /ij} 
decays. This satellite peak can be well described in MC 
simulation. In Fig. (^Id), the three Xcj tails show distin- 
guishable distributions; the small left bump is from the 
Xci tail, while the XcO tail is dominant at the right side. 
The distribution in data in Fig. (^fc) can only be well 
described by the simulated XcJ shapes. 

The angle of the normal axis of the "0(3686) decay 
plane with respect to the 0(3686) polarization vector 
(aligned to the beam axis), /3, can be determined in our 
data. The event rate may be expressed, to leading order, 
as .'^^ a oc 1 -|- a cos^ B. The measurement was carried 

a cos p ^ 

out in the rest frame of the ^^(3686) and the decay plane 
of the 0(3686) was determined with the momenta of the 
two decay particles J/tp and 7ig. The signal yields in 
each angular bin were extracted by the global fit to the 
corresponding dataset following the aforementioned pro- 
cedure. After correction of the extracted signal yields 
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TABLE II. Summary of the systematic uncertainties on the mea- 
surement of Bsig of 77 J/</) signal process, Bxc./ ^'^^ Xcj intermediate 
processes and the relative branching fractions Rmn^ following the 
notation convention in Table [J The tot systematic uncertainty is 
the square root of the sum. A dash (-) means the uncertainty is 
negligible. Values inside the parentheses are for the 'yy fi'^ fi~ mode, 
while values outside are for the 776+6" mode. Numbers without 
brackets represent uncertainties that are common to both modes. 
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FIG. 3. (a) The corrected distribution of the normal angle 
j3 of the '0(3686) decay plane, and (b) the helicity angle 9e 
of J/ip decays. The curves in (a) and (b) present the fits of 
functions Po(l + acos'^ /3) and Po{l + acos^ 6e), respectively. 



with the detection efficiency, Fig. [3ja) shows the fit to 
the distribution of | cos/3| for the sum of the two dilep- 
ton modes; we obtain a ~ 0.53 ± 0.68. 

The polarization of J/ip should be helpful in under- 
standing the mechanism of the transition process jT9| . 
The polarization parameter a can be evaluated from 
the angular distribution of the decay rate, expressed as 



dN 



cx 1 -I- acos^ 9£. Here, 



^ (with Ft and 



dcosOi^^ ' ^^^^^j " rT+2rL 

Fl being the transversely and longitudinally polarized 
decay widths, respectively) and the helicity angle 9i is 
defined as the angle of the lepton in the J/ tp rest frame 
with respect to the J/ip boost direction in the labora- 
tory frame. For fully transverse (longitudinal) polariza- 
tion, a — +1(— 1). Figure [3l[b) shows the distribution of 
I cos 6*^1 for the sum of the two dilepton modes, after cor- 
recting the signal yields for the detection efficiency and 
the lepton final state radiation effect. Our fit result is 
a = 0.08 ±0.42. 

Sources of systematic errors on the measurement of 
branching fractions are listed in Table |lll Uncertainties 
associated with the efficiency of the lepton tracking and 



identification were studied with a selected control sam- 
ple of ^^(3686) — >■ 7r+7r~(^+^~)j/^. The potential bias 
due to limiting the maximum number of photon candi- 
dates was studied by varying the limit. Throughout the 
photon energy region in this work, detection and energy 
resolution of photon are well-modeled within a 1% un- 
certainty [ill, [i^l • Detector resolution of the XcJ tails is 
taken into account up to the accuracy of the MC simula- 
tion. The corresponding systematic uncertainty is evalu- 
ated by scanning the sizes of smearing parameters within 
their errors. For the signal process, the dominant uncer- 
tainties are from the description of Xcj fine shapes, e.g., 
XcJ widths. The sensitivity to the Xcj widths is studied 
by a comparison of the signal yields based on different 
settings of the Xcj widths in modeling the Xcj resonances 
within the current world-average uncertainties. Relative 
changes of the signal detection efficiencies are assigned 
as 20%, by varying the input spin-structure within the 
measurement uncertainties and weighting the efficiencies 
in the Dalitz-like plot of Fig. [TJd). 

Many sources of systematic uncertainties in Table HIl 
cancel out when extracting the "0(3686) decay plane pa- 
rameter a and the J/ip polarization parameter a. The 
quadrature sum of the remaining systematic uncertain- 
ties are 11q'27 and 1° "il ^'^^ ^^'^ respectively. 

To summarize, the first measurement of the two- 
photon transition -0(3686) — ^ ')^J/ip was carried out at 
the BESIII experiment. The branching fraction is giv- 
en in Table HI as well as those of the cascade El tran- 
sitions. The measurement of the two-photon process 
is consistent with the upper limit obtained in Ref. @. 
The results for the signal process are presented without 
considering the possible interferences between the direct 
transition and the Xcj states, due to a lack of theoret- 
ical guidance. The distribution of the normal angle of 
the ^(3686) decay plane is characterized by the param- 
eter a = 0.53 ± 0.68(stat) j;Q'27(syst), indicating a pref- 
erence for a positive value. The J /^j polarization pa- 
rameter a was evaluated as 0.08 ± 0.42(stat)lQ*j'4(syst), 
demonstrating a competitive mixing of the longitudinal 
and transverse components. These results will help con- 
strain the strength of the coupled-channel effect in future 
theoretical calculation. The reported branching fractions 
S(V'(3686) — >■ i{1jI_1'^)xc.j) are consistent with the world 
average results [14]. The reported relative branching 
fractions of B^^ , are obtained with the world's best pre- 
cision. 
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